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ABSTRACT 

We present a high-resolution (3.6", 70 pc) 12 C0(J = 1 0) mosaic of the molecular gas in M82 
covering an area of 2.5' x 3.5' (2.8kpcx3.9kpc) obtained with the OVRO millimeter interferometer. The 
observations reveal the presence of huge amounts of molecular gas (> 70% of the total molecular mass, 
Mtot ~ 1.3 x 10 9 Mq) outside the central lkpc disk. Molecular streamers are detected in and below 
M82's disk out to distances from the center of ^1.7 kpc. Some of these streamers are well correlated 
with optical absorption features; they form the basis of some of the prominent tidal HI features around 
M82. This provides evidence that the molecular gas within M82's optical disk is disrupted by the 
interaction with M81. Molecular gas is found in M82's outflow/halo, reaching distances up to 1.2 kpc 
below the plane; CO line-splitting has been detected for the first time in the outflow. The maximum 
outflow velocity is ~ 230 km s -1 ; we derive an opening angle of ~ 55° for the molecular outflow cone. 
The total amount of gas in the outflow is > 3 x 10 8 M Q and its kinetic energy is of order 10 55 erg, about 
one percent of the estimated total mechanical energy input of M82's starburst. Our study implies that 
extreme starburst environments can move significant amounts of molecular gas in to a galaxy's halo (and 
even to the intergalactic medium). 

Subject headings: galaxies individual (M 82)-galaxies: interactions-galaxies: kinematics and 
dynamics-galaxies :ISM-galaxies: starburst— ISM: jets and outflows 



1. INTRODUCTION 

M 82 is the prototypical nearby starburst galaxy; innu- 
merous studies of M 82 have been performed in the last 
decades. Most of these studies focused on the extreme 
starburst properties and the related prominent outflow 
emerging from the inner 1 kpc of M 82 (at a distance of 
3.9 Mpc, 1' ~ 1.1 kpc; Sakai & Madore 2001). The tremen- 
dous star forming activity in M82's central region can be 
traced by the high FIR luminosity (e.g., Telesco & Harper, 
1980), the filamentary Ha outflow (e.g., Bland & Tully, 
1988; Devine & Bally, 1999), the extended X-ray outflow 
(e.g., Strickland, Ponman & Stevens, 1997) and the exis- 
tence of ~ 20 compact radio sources identified as supernova 
(SN) remnants and compact HII regions (e.g., Kronberg 
et al., 1981, Pedlar etal. 1999). The central starburst is 
believed to be fueled by the large amount of molecular gas 
present in M82's center (M tot « 2.3 x 1O 8 M , e.g., Weifi 
etal. 2001). 

This central concentration of molecular gas is very 
strong in CO emission and has been the subject of many 
studies. Most of these studies were done with single dish 
telescopes (e.g., Stark & Carlson, 1982; Young & Scoville, 
1984; Wild etal. 1992; Mao etal. 2000; Seaquist & Clark, 
2001; Taylor, Walter, & Yun, 2001) addressing both the 
excitation conditions and the large scale distribution of 



the molecular gas. Some of these studies already provided 
evidence for extended CO emission around M82's center. 
Interferometric CO studies of M 82's molecular component 
have so far been limited to the central kpc (Shen & Lo 
1993, Neininger etal. 1998, Weifi etal. 1999, Weifi etal. 
2001). In this letter, we present the first results of our ex- 
tensive mapping campaign of the central 2. 8 kpc x 3. 9 kpc 
of M82 at a resolution of 3.6" (~70pc). 

2. OBSERVATIONS 

We observed M 82 in the CO(l— >0) transition using Cal- 
tech's Owens Valley Radio Observatory's millimeter array 
(OVRO) in C, L and H configurations (10 pointing mo- 
saic). This mosaic covers a 2.5' x 3.5' field, including a 
significant fraction of the optical disk of M 82 as well as its 
prominent Ha outflow. The central pointing was centered 
near the 2.2//m peak at a = 09 ,l 55 m 52 s .2, 5 = 69°40'47.1" 
(J2000.0). The other 9 pointings were shifted by 30" 
(OVRO primary beam: 60" at 115 GHz), the field of view 
is indicated in Fig. 1 (right). Data were recorded using a 
correlator setup with 240 MHz bandwidth and 2 MHz reso- 
lution, resulting in a velocity resolution of 5.2 kms -1 with 
a total bandwidth of 600 kms -1 at 115 GHz. 

For image generation and processing we used the miriad 
software package. Because extended structure was missed 
in our interferometer data, the data cube was zero-spacing 
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corrected using the IRAM 30 m data described in WeiB 
etal. (2001). The combination was done by replacing vis- 
ibilities shorter than the shortest baseline in the final pro- 
cessed data cube by visibilities computed from the single 
dish data (this technique is described in more detail in 
Weifi etal. 2001). The final data cube presented here has 
a spatial resolution of 3.6" x 3.6" and an rms noise level 
of 25mJybeam _1 (0.18K) in a 5.2 kms -1 channel. 

3. DATA PRESENTATION AND ANALYSIS 

3.1. Global properties 

The integrated intensity distribution of our combined 
mosaic is shown in Fig. 1. The physical properties 
(CO flux, molecular mass, linear size) of the newly de- 
tected features are summarized in Tab. 1. In convert- 
ing the observed CO luminosity to H2 column densi- 
ties we applied two conversion ('X') factors 1) Am 82 = 
JV(H 2 )/I(CO) = 0.5 x 10 20 cm- 2 (Kkms^ 1 )^ 1 , an aver- 
age value for M82's central disk (Weifi etal. 2001), and 2) 
^Gai = 1.6xl0 20 cm~ 2 (Kkms -1 ) -1 , the Galactic conver- 
sion factor (Hunter et al. 1997) which is presumably a more 
realistic value for a non-starburst environment (such as 
the molecular streamers S1-S4, see Fig. 1, right). The to- 
tal flux detected in the field of view is 17680 Jy km s -1 (es- 
timated error: ~ 10%). Using Am 82 for the disk and the 
outflow O-S + O-N (which presumably consists of warm 
gas emerging from the disk, see Fig. 1, right) and Acai for 
the molecular streamers S1-S4, we derive a total molecular 
mass of M82 of ~ 1.3 x 10 9 M Q (total masses using other 
conversion factors can be easily derived from Tab. 1). 

3.2. Seperating the CO Components in M 82 

The different CO emitting components in M82 (disk, 
streamers, outflow/halo) have clear kinematic signatures 
which can be seperated by inspecting consecutive chan- 
nel maps. We used this velocity information to subdivide 
the full OVRO CO data cube into three sub-cubes out of 
which individual integrated intensity maps were created 
for the disk, streamer and outflow/halo component. These 
components are shown in Fig. 3 (overlaid as countours on 
an Ha map by Shopbell et al. 1998): The yellow contours 
represent M82's central disk only (see Sec. 3.3); streamers 
S1-S4 are shown in orange (Sec. 3.4), and the red contours 
represent an intensity map of the outflow/halo (Sec. 3.5). 

3.3. The Central Molecular Disk 

The central 1 kpc of M 82's CO distribution is domi- 
nated by the well known triple-peaked molecular distri- 
bution (Fig. 1, see also yellow contours in Fig. 3) which 
has been discussed in previous high-resolution CO studies. 
This region coincides with the strong central star forma- 
tion and forms the base of the prominent outflow visible 
in Ha and X-rays. It is striking that only about a third 
of M82's total CO flux is associated with the center (see 
Tab. 1). Depending on the employed X-factor this cor- 
responds to a molecular mass fraction of < 30% for the 
central molecular disk. 

3.4. Molecular Streamers 

Several molecular streamers are detected west, east and 
south-east of the central molecular concentration (labeled 



SI to S4 in Fig. 1, see also orange contours in Fig. 3). The 
most prominent streamer is located west of the molecular 
disk (SI). Its orientation is shaped in a similar fashion and 
has the same velocity as the prominent H I streamer which 
points towards M81 (the M81-M82 tidal tail, streamer 'S' 
in Yun, Ho & Lo 1993). A weaker streamer forms the 
extension of the molecular disk towards the north-east 
(S2), roughly following another large-scale HI streamer 
(streamer 'N' in Yun, Ho & Lo 1993). SI and S2 follow 
very closely the dust absorption features seen in optical 
images (Fig. 2). Two additional molecular streamers are 
visible south-east (S3) and east (S4); in contrast to SI 
and S2, these two structures are situated away from the 
optical plane of the galaxy (see also Figs. 1, 3 and 5, bot- 
tom). Based on the available data, we can not rule out a 
potential interaction between S3/S4 and the outflow/halo 
gas. A summary of the fluxes of S1-S4 is given in Tab. 1. 



3.5. Molecular Gas in the Outflow/ Halo 

Masses: Molecular gas is clearly associated with the out- 
flow south (O-S) and north (O-N) of the central molecular 
disk. The molecular gas in the outflow/halo is visible up 
to 0.8 kpc north and 1.2 kpc south of the disk. The total 
mass of molecular gas in the outflow/halo accounts to an 
impressive amount of 3.3 x 10 8 M Q , even employing the 
lower conversion factor Xm82 ■ 

Outflow morphology: Interestingly, the CO emission in 
the outflow/halo (Fig. 3, red contours) has an almost 
spherical shape with a diameter of <~ 1.5 kpc (this is simi- 
lar in size to the dust halo discussed in Alton et al. 1999). 
In addition, there is a striking correlation between the Ha 
outflow and the molecular gas in the south. The corre- 
lation is less prominent in the north-east of the outflow 
but clearly visible in the north-west. The chimney of 
shocked molecular gas detected in SiO by Garcia-Burillo 
et al. (2001) can not be clearly identified in our CO data. 

Kinematics: For the first time CO line splitting is de- 
tected in the outflow of M 82. Fig. 4 shows a pV diagram of 
the CO gas along the outflow axis of M 82. Here, a position 
angle of 150° was chosen to allow a comparison with pub- 
lished Ha velocities (see caption). The CO line splitting 
is clearly visible between —20" and —35" below the disk; 
the velocities agree well with the Ha velocities. The blue- 
shifted CO component is not detected below —35" but the 
other component is traced at similar velocities as the Ha 
up to -70". At -35" (660 pc below the plane), the CO is 
blucshifted by <~ 130kms~ 1 with respect to the systemic 
velocity (v sys =220kms~ 1 ). Correcting for M82's inclina- 
tion (80°, de Vaucouleurs et al. 1991) and our derived cone 
angle of 50° (see below) we derive a deprojected maximum 
outflow velocity of ~ 130kms _1 /sin(35°) w 230 kms -1 . 
The situation is more complex in the northern part of the 
outflow - here the molecular gas seems to be entrained 
in the Ha outflow. As is obvious from Fig. 4, most of 
the projected velocities are lower than this value and we 
estimate an average outflow velocity of the gas of order 
~ 100 kms -1 . Based on the escape velocities presented in 
Martin (1998, Figure 9) we estimate that of order 10% of 
the outflow gas (3.3 x 10 7 M Q , i.e., a few percent of the 
total gas mass) may be lost to the intergalactic medium 
(IGM). 
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Outflow cone angle: Fig. 5 presents pV diagrams along 
and parallel to the major axis. In the north, the rotation 
curve is flattened which can be interpreted as conserva- 
tion of angular momentum in the outflow. A similar mor- 
phology is seen towards the southern pV diagram as well; 
however, this region is clearly disturbed by the features S3 
and S4. 

To determine the orientation of the outflow cone, we 
derived 'flow-lines' of the molecular gas (as first proposed 
by Seaquist & Clark 2001). The basic assumption is that 
the angular momentum of the gas in the outflow is con- 
served. For each distance from the plane z we calculated 
the angular momentum L ~ r ■ v t where ft is the termi- 
nal velocity at distance r from the minor axis. We took 
Vt from our first moment map of our 'outflow cube' (see 
above); these values give a reasonable approximation for 
the true v t (see also the discussion in Seaquist & Clark 
2001). The result of this analysis is plotted as contours in 
Fig. 6 - each contour represents an iso-angular momen- 
tum ('-£<') line and describes the flow-line of an individual 
particle in the outflow. It is striking how similar the L- 
contours (derived entirely from kinematic information) are 
compared to the morphology of the ionized gas in the out- 
flow (e.g., Fig. 3). From the third contour (counting from 
the minor axis) in Fig. 6, we derive an average cone angle 
of ~ 55° (north: ~ 60°, south: ~ 50°), similar to what has 
been found in previous optical studies (optical: <~ 30 — 35°, 
Bland & Tully 1988, McKeith etal. 1995, to - 60 - 65°, 
Heckman, Armus & Miley 1990). We note that Seaquist 
& Clark (2001) derived an angle of ~ 40° for the northern 
and ~ 90° for the southern molecular outflow; the discrep- 
ancy in the south can be explained by confusion of their 
velocity field by S3 and S4 in their low-resolution data. 

Energetics: To calculate an order of magnitude estimate 
of the kinetic energy of the outflow, we assume that the 
molecular gas in the outflow (3.3 x 1O 8 M ) has an average 
outflow velocity of ~ 100 kms -1 (see discussion above). 
This gives a kinetic energy of ~ 3.3 x 10 55 erg. Extrapolat- 
ing the current SN rate in M 82 of 0.1 yr _1 (Kronberg et al. 



1981) back to the age of the current starburst (5 x 10 7 yr, 
O'Connell etal. 1995) yields a total mechanical input by 
SNe of <~ 5 x 10 57 erg. Assuming an average outflow veloc- 
ity of 100 km s _1 and a distance of ~ 1 kpc it takes of order 
1 x 10 7 yr to build up the molecular outflow, similar to the 
age of the starburst. This implies that of order 1% of the 
released mechanical energy of M82's central starburst has 
been converted to kinetic energy in the molecular outflow 
gas. 

4. CONCLUDING REMARKS 

Our high-resolution OVRO observations of the proto- 
typical starburst galaxy M 82 provide new exciting insights 
on the properties of the molecular gas around the galaxy's 
center. The newly detected molecular streamers (both in 
and below the galaxy's plane) suggest that the molecu- 
lar gas in the inner few kpc is severely affected by the 
interaction with M81. This redistribution of molecular 
gas is likely the trigger for the strong starburst activity in 
M82's center. An impressive amount of molecular gas is 
detected in the halo of M 82. It is striking that most of 
the outflow gas is distributed in an almost spherical halo; 
at large distances from the plane its morphology looks re- 
markably similar to that of the ionized gas. The gas in 
the halo is rotating slower than in the disk; we find that 
the CO kinematics in the halo can be used to constrain 
the outflow cone angle of the molecular gas (<~ 55°). Our 
estimated CO outflow velocities imply that a significant 
amount of molecular gas is not only transported to the 
halo but may also escape the galactic potential. This has 
important implications for the chemical evolution of the 
starburst host galaxy and for the chemical enrichment of 
the IGM, especially in the early universe when much more 
violent starbursts took place. 

FW acknowledges NSF grant AST 96-13717. AW ac- 
knowledges DFG grant SFB 494. We thank P. Shopbell 
for providing us with his M 82 Fabry-Perot Ha data. 
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Table 1 
Observed CO fluxes in M82 
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Fig. 1. — Left: Logarithmic representation of the integrated CO(l— >0) map of the zero-spacing corrected OVRO mosaic (moment 0). 
Contours arc plotted at 1.5, 3, 6, 12, 24, 48, 96, 192 Jy beam" 1 krns" 1 ; the rms is 0.3 Jy beam 'kms 1 . The OVRO mosaicked field of 
view is indicated by the outer envelope (plotted at 1.5 X the FWHM of the primary beam). Right: Same plot with superimposed labels of 
the molecular streamers (S1-S4) and the outflow gas (O-N and O-S); the lines indicate the orientations of the position velocity diagrams 
discussed in Fig. 4 and 6. The dashed lines indicate the derived opening angle of the molecular outflow (50°, south; 60°, north, see Sec. 3.5). 



Fig. 2. — Three color composite of M82: Red: Optical image of M82; Blue: Ha image of M82 (see Fig. 3); Green: OVRO CO map (see 
Fig. 1). 



Fig. 3. — H« image of M82 (greyscale, see text). White contours: M82 disk, Grey contours: M82 streamer (S1-S4), Black contours: 
molecular gas associated with the outflow (O-S + O-N). Contour levels are the same as in Fig. 1, a contour of 0.4 Jy beam -1 kms -1 has been 
added to the outflow component. Note that molecular gas is clearly associated with M82's prominent outflow of ionized gas. 



Fig. 4. — CO pV diagram along the outflow in M82 (PA=150°, see 'pV-O' in Fig. 1 for orientation). The crosses (Shopbell etal. 1998) and 
pluses (McKeith etal. 1995) represent the Ho velocities. Towards the south (negative offsets), CO line splitting is detected. 



Fig. 5. — CO pV diagram along and parallel to the major axis of M82 (PA=75°, see labels 'pV-C, 'pV-N' and 'pV-S' in Fig. 1 for the 
relative orientations, offsets are given with respect to the minor axis). Contours in the south arc shown at 0.08, 0.15, 0.2, 0.25 Jybeam - 1 , 
along the major axis contours are also plotted at 0.75, 1.5Jybeam . Note that the rotational amplitude is decreasing with distance from 
M82's plane. The southern diagram (bottom) includes part of the streamers SI, S3 and S4. 



Fig. 6. — Iso-angular momentum ('!/') contours (black) of molecular gas in M82 (greyscale: moment map, same as in Fig. 1). Under the 
assumption that angular momentum in the outflow is conserved, these contours show the 'flow-lines' of the molecular gas (see text). 
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